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SUMMARY 

^ II ') 

/  ^  f-By  the  use  of  orthogonal  polynomials  developed  for  discrete  sets  of  data, 

/  the  least- squares  equations  for  determining  the  optimized  stress-rupture  para- 
I  metric  constants  are  obtained  in  nearly  uncoupled  form;  thus  the  use  of  high- 

I  degree  polynominals  is  permitted  without  the  loss  of  significant  figures. 

/  Optimum  values  of  the  constants  can  thereby  be  accurately  obtained.  The  method 

'  is  applied  to  the  data  obtained  from  the  German  cooperative  long-time  creep 

program  by  using  a  general  parameter  of  which  the  Manson-Haferd  and  Larson- 
Miller  parameters  are  special  cases.  Good  correlation  was  obtained.  An  analy¬ 
sis  is  also  made  of  creep  data  obtained  for  columblum  alloy  FS-85  with  good 
results.  A  complete  Fortran  IV  computer  program  is  included  to  aid  those  wish¬ 
ing  to  use  the  method.  I 

I  -^^s" 
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INTRODUCTION 
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One  method  of  extrapolating  short-time  creep-rupture  data  to  predict  long¬ 
time  life  involves  the  use  of  a  time-temperature  parameter.  This  concept  is 
based  on  the  assumption  that  all  creep-rupture  data  for  a  given  material  can  be 
correlated  to  produce  a  single  "master  curve"  wherein  the  stress  (or  log  stress) 
is  plotted  against  a  parameter  involving  a  combination  of  time  and  temperature. 
Extrapolation  to  long  times  can  then  be  obtained  from  this  master  curve,  which 
can  presumably  be  constructed  by  using  only  short-time  data.  Three  well-known 
parametric  methods  are  the  Larson-Miller,  Manson-Haferd,  and  Dorn  parameters 
(refs.  1  to  3).  These  parametric  methods  have  the  great  advantage,  at  least  in 
theory,  of  requiring  only  a  relatively  small  amount  of  data  to  establish  the 
required  master  curve. 

All 

More  recently  a  general  creep- rupture  parameter  was  introduced  by  one  of 
the  authors  (ref.  4)  that  includes  most  of  the  currently  used  parameters  as 
special  cases.  The  analysis  in  the  present  paper  is  therefore  based  on  this 
general  parameter. 
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A  significant  advance  in  the  practical  application  of  the  parametric 
methods  was  the  development  of  an  objective  least-squares  method  for  determin¬ 
ing  the  optimum  values  of  the  parametric  constants  without  plotting  and  cross - 
plotting  the  data  and  without  the  use  of  judgment  on  the  part  of  the  analyst 
(ref.  5).  This  least-squares  method  involves however^  several  practical  dif¬ 
ficulties  that  arise  from  the  fact  that  in  fitting  the  master  curve  by  a  poly¬ 
nomial^  the  set  of  linear  algebraic  equations  for  the  coefficients  (the  normal 
equations)  are  very  ill-conditioned.  The  determinant  of  these  equations  can  be 
shown  to  be  related  to  the  Hilbert  determinant  (ref.  6)^  which  rapidly  ap¬ 
proaches  zero  as  its  order  increases.  Thus  for  polynomials  above  the  second 
degree;  it  is  necessary  to  use  double-precision  arithmetic  (l6  significant 
digits  or  more)  on  the  computer;  and  for  the  fifth  degree  and  above  the  results 
become  uncertain  even  with  double-precision  arithmetic.  This  difficulty  is 
inherent  in  the  normal  least-squares  equations  and  is  not  limited  only  to  the 
stress-rupture  problem. 

The  present  report  presents  a  method  for  avoiding  the  above  difficulty  by 
using  orthogonal  polynomials  in  the  representation  of  the  master  curve  (appen¬ 
dix  A).  The  use  of  orthogonal  polynomials  for  representing  discrete  sets  of 
unequally  spaced  data  is  described  in  reference  6  and  in  more  detail  in  refer¬ 
ence  7.  A  further  improvement  can  be  obtained  by  performing  a  linear  trans¬ 
formation  on  the  stresses  (or  the  logs  of  the  stresses)  so  that  all  the  values 
of  stress  (or  log  stress)  lie  between  2  and  -2;  as  recommended  in  reference  7. 

As  a  result  of  these  innovations;  it  became  possible  to  perform  all  the  compu¬ 
tations  in  single-precision  arithmetic  (eight  significant  digits)  up  to  18th 
degree  polynomials  without  appreciable  round-off  error. 

In  addition;  this  report  contains  a  complete  analysis;  in  which  the  gen¬ 
eral  parameter  was  used;  of  all  the  data  for  three  steels  that  were  obtained  by 
NASA  through  the  cooperation  of  Dr.  K.  Richard  of  Faberwerke  Hoechst  in 
Frankfurt  and  that  were  investigated  in  a  long-time  cooperative  creep  program 
in  Germany.  Some  of  the  data  from  the  latter  investigation  are  included  in  this 
paper. 

Finally  it  is  shown  by  means  of  a  concrete  example  how  the  parameter  tech¬ 
niques  can  be  applied  to  creep  data  to  predict  long-time  creep.  For  this  pur¬ 
pose  the  data  for  columbium  alloy  FS-85;  as  reported  in  reference  8;  are  used. 

A  complete  Fortran  IV  program;  as  used  on  the  IBM  7094  computer  in  making 
the  calculations;  is  presented  in  appendix  B.  This  program  can  be  used  for  the 
objective  analysis  of  any  set  of  creep-rupture  data  by  the  Larson-Miller; 
Manson-Haferd;  or  the  more  general  parameter  of  reference  4. 


SYMBOLS 

A;B  linear  transformation  coefficients 

a^b;C  elements  of  coefficient  matrix 

D  standard  deviation 
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K  degree  of  freedom 

m  degree  of  polynomial 

n  nnmler  of  data  points 

P(a)  creep-rupture  parameter 
Q  polynomial 

q  stress  exponent 

r  temperature  exponent 

S  sum  of  squares  of  residuals 

T  temperature 

temperature  intercept 
t  time  to  rupture 

tg,  time  intercept 

u  coefficient  of  polynomial  function 

X  scaled  log  stress 

X  log  stress 

y  log  time 

yg^  log  time  intercept 

constants  from  recurrence  relation 
0  stress 

T  a'lCT  -  Ta)^ 

Subscripts : 
max  maximum 

min  minimum 

PROCEDURE 
General  Parameter 

The  general  creep-rupture  parameter  introduced  in  reference  4  has  the  fol- 
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lowing  form 


P(a) 


-  log  t 

_o^ _ ^ 

(T  -  Ta)^ 


(1) 


where  log  tg,  q,  and  r  are  material  constants  to  he  determined  from  the 

available  experimental  data.  The  parameter  P(ct)  is  a  function  of  the  stress 
and^  when  plotted  against  stress,  is  referred  to  as  a  master  curve  (fig.  1, 
p.  9) .  If  q  =  0  and  r  =  1,  the  Manson-Haferd  parameter  is  obtained.  If 

q  =  0,  r  =  -1,  and  =  -460^  F,  the  Larson-Miller  parameter  results.  If 

q  =  1  and  r  =  1,  the  stress-modified  parameter  suggested  in  reference  9  is  ob¬ 

tained.  Finally,  if  q  =  0,  equation  (l)  reduces  to  the  parameter  proposed  by 
Manson  and  Brown  (ref.  10). 


The  object  is  to  find,  the  best  values  of  the  constants  q,  log  t^^  Tg^, 
and  r  so  that  the  master  curve  best  fits  the  data.  To  find  these  values,  the 
method  of  least  squares  is  used  whereby  the  master  curve  is  represented  by  a 
polynomial  in  the  logarithm  of  the  stress,  and  the  best  fit  is  obtained  by  mini¬ 
mizing  the  sum  of  the  squares  of  the  deviations  (the  residuals)  of  the  data 
from  the  curve.  The  calculation  procedure  will  now  be  described.  The  details 
of  the  derivation  are  given  in  appendix  A,  and  a  Fortran  IV  computer  program 
using  this  method  is  given  in  appendix  B. 


Calculation  Procedure 


To  simplify  the  notation,  the  following  symbols  are  introduced: 

T  =  ct'1(T  -  Ta)^ 
y  s  log  t 

>• 

X  =  log  a 
Ya  =  log  la 


(2) 


Then  from  equation  (l)  it  follows  that 

y  =  o'S'a  +  'rQ(x)  (3) 

where  in  reference  5,  Q(x)  was  represented  hy  a  simple  polynomial  of  the  form 


Q(x)  =  bq  +  a^x  +  agx^  +  .  .  .  +  a^jX"^  (4) 

The  least-squares  equations  obtained  sometimes  led  to  difficulties  as  indicated 
in  the  lOTRODUCTION.  These  difficulties  can  be  avoided,  however,  by  rewriting 
equation  (4)  in  terms  of  polynomials  that  are  orthogonal  over  the  set  of  data, 
as  defined  in  appendix  A.  Thus  assume 
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Q(x)  =  Uj^Q^(x)  +  UgQ^Cx)  +  .  . 


(5) 


m+1 

•  +  “m+A+l'=‘>  “  E 

j=l 

■where  u-  is  an  unkno-wn  constant^  m  is  the  degree  of  the  highest  degree 
polynomial^  and  Qj(x)  is  a  polynomial  of  degree  j  -  1  that  satisfies  the 
orthogonality  conditions  described  in  appendix  A.  The  use  of  orthogonal  poly¬ 
nomials  permits  the  solution  of  the  least-squares  equations  directly  in  closed 
form,  thus  the  loss  of  a  large  number  of  significant  digits  is  avoided.  The 
method  of  calculating  Qj  ’will  be  discussed  in  appendix  A. 

If  equation  (s)  is  substituted  into  equation  (3),  an  equation  -with  m  +  5 
unkno-wn  constants  results  for  the  case  of  the  general  parameter.  For  the  case 
of  the  linear  parameter  there  are  m  +  3  constants ,  and  for  the  Larson-Miller 
parameter  there  are  m  +  2.  It  is  necessary  that  the  number  of  data  points  n 
always  equals  or  exceeds  the  nurtiber  of  unknown  constants. 

The  constants  are  deteimined  so  that  equation  (3)  fits  the  data  best  in 
the  least-squares  sense.  To  accomplish  this,  the  sum  of  the  squares  of  the 
deviations  is  minimized;  that  is, 

S  ^  E  &i  ■  ^i^a  -  '^i^(^iO^ 

1=1 

is  made  a  minimum.  Because  the  equations  are  nonlinear  in  some  of  the  unkno'wn 
constants  a  trial  and  error  procedure  must  be  used.  A  set  of  values  is  assumed 
for  q,  r,  and  Ta^  and  the  corresponding  best  values  of  ya  a,nd  uj  are  de¬ 
termined.  A  different  set  of  values  for  q,  r,  and  Tg,  is  then  chosen,  and 
again  the  best  values  of  yg  and  Uj  are  calculated.  Several  sets  of  -values 
of  q,  r,  and  Tg  are  tried,  and  the  values  corresponding  to  the  overall  best 
fit  are  determined.  For  the  case  of  the  linear  parameter,  only  the  value  of 
Ta  is  varied  (q  is  always  equal  to  zero,  and  r  is  always  equal  to  l).  For 
the  Larson-Miller  parameter,  Tg  is  equal  to  -460°  F,  and  no  trial  and  error 
procediire  is  needed. 

As  a  measure  of  the  fit,  the  standard  de-viation  D,  defined  by 
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To  determine  the  test  values  of  and  uj  for  a  given  set  of  values  of 

following  calculations  are  made.  Firsts  the  logarithms  of 
the  stresses  are  scaled  so  that  they  lie  in  the  range  -2  to  2^  as  suggested  in 
reference  7.  The  reason  for  this  is  discussed  in  appendix  A.  Thus  define  a 
variable  X  by 


A  =: 


B  =  -2 


X  =  Ax  +  B 
4 

^ax  “  ^in 
^max  %in 


> 


^ax  “  ^in 


(9a) 


(91d) 


J 


The  polynomials  Qj(Xj_)  are  now  calculated  for  each  of  the  data  points  by 
using  the  following  formulas: 


m  >  j  >  1 


(10) 


n 


2  XiT?Q2(Xi) 


i=l 


a.  =  - 

J  n 


m  >  j  >  1 


X; 


i=l 


n 


y  (lOa) 


i=l 


n 


m  >  j  >  1^  Q-,  =1;,  and  P-.  =  0 


2  ^iQj-i(Xi) 


i=l 


where  n  is  the  number  of  data  points^  X^^  is  the  scaled  value  of  log  for  the 
i"^^  data  pointy  and  Tj_  is  equal  to  data  point  for 

the  chosen  values  of  Tg^^  q^,  and  r. 

It  is  to  be  noted  that  the  degree  of  the  polynomial  Q(x)  of  equation  (5) 
can  be  increased  by  merely  computing  the  next  polynomial  in  the  series  *^.^.2 

without  having  to  recompute  any  of  the  previous  ones.  This  is  one  of  the  ad¬ 
vantages  of  using  orthogonal  polynomials. 

Once  the  values  of  Qj  have  been  computed  for  each  of  the  data  points^ 
yg^  and  Uj  can  be  calculated  as  follows: 
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Let 


■where  j  =  1^2 
Then 


n 


QO 


2q 


i=l 


n 


i=l 


n 


”0  ■  E 


i=l 


n 


'0 


=  E 

1=1 


^lyi 


n 


j  =  X) 


i=l 


m  +  1. 


= 


m+1 

'0  Aj  b. 
j=l  J 


m+1 


'0 


.  y  f:i 


j=i 


y  (11) 


>  (12) 


u 


J 


■ 


Note  that  if  q.  =  0^  equals  the  number  of  data  points  n.  Thus  by  means  of 
equations  (9)  to  (12)>  the  test  values  of  ya  and  uj  to  fit  the  data  are 
found  for  a  given  choice  of  Tg^^  q^  and  r.  The  Fortran  IV  program  described 
in  appendix  B  automatically  scans  all  the  desired  values  of  Ta^  q^  and  r  and 
chooses  the  best  set  from  all  the  submitted  values  as  determined  by  the  smal¬ 
lest  value  of  the  standard  deviation  as  defined  by  equation  (?).  The 
method  can  be  illustrated  by  a  simple  example:  consider  a  set  of  theoretical 
data^  which  fit  the  following  equation  exactly 
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(13) 


Q-5  -  ,1.9.g_t  =  10“^(7.02  +  0.467  x  +  0.061  x^  +  0.00928  x^) 

T  -  600 

Eight  data  points  satisfying  this  equation  are  given  in  colianns  2  to  6  of 
table  I.  For  this  data  Ta  =  600°  F  and  log  ta  =  Ya  =  9* 6*  Suppose^  however^ 
that  these  eight  data  points  were  obtained  experimentally  and  that  the  values 
of  Ta  s-nd  log  tg^  were  not  known.  The  problem  then  is  to  find  the  best 
values  of  Ta  and  log  ta  to  fit  the  data  by  the  linear  parameter.  These 
values  can  readily  be  found  by  using  the  equations  of  the  previous  section. 
Firsts  from  column  6  of  table  I 

(log  =  4.  75051 

(log  =  1.81954 


Therefore  from  equations  (9b) 


A  =  1. 56474 
B  =  -4.48319 

and  by  means  of  equation  (9a)  the  were  computed  and  are  given  in  column  8. 

For  illustrative  purposes  three  values  of  Ta  were  chosen,  500°,  600°, 
and  700°  F.  For  each  of  these  values  of  Ta,  values  of  Ti,  a^,  pj,  and 
Sj(Xl)  were  computed  by  means  of  equations  (2),  (lO),  and  (lOa;,  and  the  values 
of  aj,  bj,  and  Cj  were  computed  by  equations  (ll).  The  results  are  tabu¬ 
lated  for  Ta  =  600°  in  columns  9  to  12  of  table  I  and  in  table  II  up  to  a 
third  degree  polynomial. 

The  values  of  ya  and  Uj  were  then  computed  by  using  equations  (l2)  for 
each  of  these  three  values  of  Ta  by  first  assuming  m  =  2,  then  m  =  3,  and 
finally  m  =  4,  corresponding  to  polynomials  of  second,  third,  and  fourth  de¬ 
grees,  respectively.  For  each  of  these  cases  the  standard  deviation  D  was 
computed  from  equation  (7)  with  S  being  given  by  equation  (6)  and  Q  by 
equation  (5).  The  results  are  simmarized  in  table  III.  The  least  value  of  D, 
signifying  the  best  fit,  is  obtained  for  m  =  3  and  Ta  =  600°  F.  The  cor¬ 
responding  value  of  •  yg  is  9.5.  These  values,  of  course,  correspond  to  equa¬ 
tion  (13),  from  which  the  data  were  generated. 


Application  to  Data  from  German  Cooperative  Long-Time  Creep  Program 

As  part  of  the  German  cooperative  long-time  creep  program,  a  sufficient 
amount  of  material  of  each  of  three  steels  was  supplied  to  NASA  to  permit  the 
running  of  short-time  tests  necessary  to  predict  the  results  at  long  times  ob¬ 
tained  in  the  German  test  program.  The  composition  of  these  steels  is  shown  in 
table  IV. 

The  results  of  the  NASA  tests,  which  were  used  in  the  subsequent  analysis. 
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are  shovn  in  table  V.  Table  VI  shows  the 
results  of  the  long-time  German  test  pro¬ 
gram.  The  three  steels  will  be  designated 
briefly  as  steel  steel  and  steel  P. 

With  the  use  of  the  test  data  shown 
in  table  V  a  complete  analysis  was  made 
by  the  previously  described  method.  The 
general  parameter  discussed  in  the 
IWTRODUCTIOW  was  used^  arid  the  best  values 
were  obtained  for  the  parametric  constants 
for  each  of  the  three  steels. 


All  the  data  obtained  for  these 
steels  are  shown  in  tables  V  and  VI. 

Many  of  the  data  points  were  obtained  for 
purposes  other  than  the  application  to 
time-temperature  parameters^  as  described 
in  this  report.  As  already  discussed  in 
references  4  and  11^  a  much  smaller 
amount  of  data  is  needed  when  an  accel¬ 
erated  program  is  desired;  however^  since 
these  data  were  already  available^  all  the 
data  indicated  in  tables  V  and  VI  were 
used  to  obtain  the  best  possible  para¬ 
metric  constants. 

For  all  three  steels  the  analysis  showed  the  stress  exponent  to  be 
zero^  but  the  temperature  exponent  r  to  be  different  for  each  of  the  three 
materials.  For  steel  K  the  best  value  of  r  was  1^  which  indicated  that  the 
best  fit  is  obtained  by  the  linear  parameter.  For  steel  P  a  value  of  r  of 
-1  was  obtained^  which  indicated  a  parameter  similar  to  the  Larson-Miller  param¬ 
eter;  however^  the  corresponding  value  of  Ta  tos  200^  F  rather  than  -460^  F 
used  in  the  Larson-Miller  parameter.  For  steel  C  the  value  of  R  was  2.5. 

Figure  1  shows  the  results  for  steel  K.  Here  the  master  curve  consists  of 
a  plot  of  stress  against  the  optimized  parameter  (T  -  300)/(log  t  -  16.54). 

Figure  2  shows  the  isothermals  computed  by  using  the  optimized  parameters^ 
as  shown  on  each  of  the  figures.  The  range  of  the  NASA  data  used  to  obtain 
these  parameters  is  also  shown  on  each  of  the  figures.  The  data  points  shown 
are  the  German  results  obtained  to  date.  G?he  predictions  up  to  100  000  hours 
from  the  NASA  data  based  on  the  optimized  parameters  agree  well  with  the 
German  data,  if  scatter  and  differences  in  testing  technique  between  the  two 
organizations  are  considered. 

Figure  3  shows  a  comparison  for  each  of  the  three  steels  between  the  best 
linear  parameter,  the  best  Larson-Miller  parameter,  and  the  best  general  param¬ 
eter.  Although  for  some  of  the  steels  fair  agreement  can  be  obtained  with  one 
or  the  other  of  these  parameters,  it  is  clear  that  the  general  parameter  is 
superior  when  all  the  materials  are  considered  jointly.  If  any  one  of  the 
special  cases  of  this  parameter  is  to  be  chosen  for  all  materials,  the  linear 


Figure  1.  -  Master  curve  for  steel  K  (27b  KK),  calculated 
from  NASA  data  between  10  and  3700  hours. 
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Stress,  ksi 


(a)  Steel  K  (27b  KK);  parameter,  (b)  Steel  C  (23b  CK);  parameter, 

P  =  (log  t  ‘  16. 54)/(T  -  300).  P  =  (log  t  -  8. 87)/(T  +  100)  •  ^ 

Figure  2.  -  Analysis  of  German  steel  data  by  generalized  parameter  with  optimuip  constants  (where 


(c)  Steel  P(14a  PA);  parameter, 

P  =  (T-  200)(logt  +  11.13). 

T  is  temperature,  and  t  is  time  to  rupture). 


Temperature  for 

German  data 
above  5  hr. 

°F 

Q 

932 

o 

1022 

o 

1112 

A 

1202 

□ 

1292 

Figure  3.  -  Analysis  of  German  steel  data  by  several  parameters  (where  T  is  temperature,  and  t  is  time  to  rupture). 
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Figure  4.  -  Analysis  of  creep  data  for  columbium  alloy  FS“85  by  linear  parameter. 


parameter  vould  appear  to  be  the  best  choice. 


Application  to  Creep  Data 

^^though  there  is^P,  fundamental  reason  why  the  same  parameter  is  capable 
of  representing  both  creep  and  rupture  data^  it  has  nevertheless  been  found 
empirically  (refs.  1  and  2)  that  the  dual  role  of  the  same  parameter  leads  to 
reasonable  results.  Experimental  data  for  creep  are  much  more  limited,  however, 
than  that  for  rupture,  and  such  data  tend  to  contain  more  scatte^  hence, 
analysis  of  creep  data  by  the  parametric  approach  has  been  limited'  in  the  past. 

The  method  of  the  present  report  can  be  applied  directly  to  creep  data 
without  any  change.  All  that  is  necessary  is  to  redefine  t  as  the  time  to 
attain  a  specified  amount  of  creep  rather  than  as  the  rupture  time.  Thus,  it 
is  assumed  that  for  a  given  amount  of  creep,  say  1  percent,  a  plot  of  log  o 
against  a  parameter,  such  as  that  given  by  equation  (l) ,  will  produce  a  single 
master  curve.  For  a  different  amount  of  creep,  say  5  percent,  a  different 
master  curve  can  be  obtained,  but  it  is  assumed  that  the  parametric  constants, 
such  as  log  ta  and  Tg^,  remain  the  same  and  that  they  equal  the  values  ob¬ 
tained  from  rupture  data. 

Calculations  of  this  type  were  performed  for  columbium  alloy  FS-85.  The 
creep  tests  were  limited  to  runs  of  approximately  1000  hours,'  the  data  are 
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given  in  table  'VII,  as  taken  from  reference  8.  Figure  4(a)  shows  the  data  for 
5-percent  creep  strain,  and  figure  4(b)  shows  the  master  curves  obtained  for 
1-,  2-,  and  5-percent  strain  as  well  as  the  parametric  constants  obtained  by 
the  method  of  this  report.  While  scatter  in  the  creep  data  is  high,  the  cor¬ 
relation  must  be  regarded  as  good.  In  general,  the  points  agree  well  with  the 
master  curve. 

Although  these  results  are  encouraging,  much  more  work  is  necessary  before 
it  can  be  concluded  that  the  parametric  approach  is  completely  valid  for  creep 
data.  If  it  is  eventually  concluded  that  the  parametric  approach  is  valid  for 
creep  data  and  in  particular  that  the  parametric  constants  are  the  same  for 
both  the  creep  and  rupture  processes,  it  is  obvious  that  a  great  saving  in  test 
facilities  and  test  program  planning  will  result.  It  therefore  seems  very 
worthwhile  in  future  studies  to  give  more  attention  to  the  correlation  and 
extrapolation  of  creep  data  by  the  parametric  method. 


Lewis  Research  Center, 

National  Aeronautics  and  Space  Administration, 
Cleveland,  Ohio,  May  3,  1965. 
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APPEIDIX  A 


ORTHOGOKAL  POLYNOMIALS  AND  LEAST- SQUARES  DETERMINAriON 
OP  PARAMETRIC  CONSTANTS 

A  set  of  polynomials  Qj(x)  are  said  to  "be  orthogonal  over  an  interval 
with  respect  to  the  weighting  function  t(x)  if  they  satisfy  the  following 
relation 


j.  k  (Al) 


Similarly  a  set  of  polynomial’s  can  be  defined  to  be  orthogonal  over  a  set  of  n 
discrete  points  by  the  following  relation 

n 

2  TiQj(xi)Qjj(Xi)  =0  j  j/  k  (A2) 

i=l 

It  can  be  shown  (ref.  6)^  that  all  orthogonal  polynomials  satisfy  a  three-term 
recurrence  relation  of  the  form 

^k+1  =  (^  ‘  "  Pk*^k-1  ^  ^  ^ 

Thus  by  starting  with  =  1  and  =  0  an  infinite  set  of  orthogonal  poly¬ 
nomials  can  be  generated  by  means  of  equation  (A3)  if  values  for  and 

are  known.  These  can  be  determined  from  the  orthogonality  conditions  (eqs. 

(Al)  or  (A2)).  From  the  relation  (A2)  it  follows  that 

n 

^T?Qk(Xi)Qk+i(xi)  =  0  (A4a) 


n 

^  '^i\+l(^i)\-l(^i)  =  ° 

i=l 

When  the  recurrence  relation  (A3)  is  used  to  eliminate  there  is  obtained 

n 

D  -  f>A-i]  -  0 

i=l 
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i=l 


(ASTo) 


■When  the  orthogonality  condition  (A2)  is  used^,  equations  (A5a)  and  (ASh)  reduce 
to 


n 


s 

1=1 


|(xi  -  a^)Q|  =  0 


(A6a) 


n 


1=1 


Solving  equations  (A6)  for  and  gives 


n 


E2„2 

XiTiQj, 


1=1 
°^k  =  — 


2n2 


1=1 


(A6Td) 


(A7a) 


n 


i=l 


n 


E 


1=1 


(A7b) 


Thus  a  set  of  orthogonal  polynomials  can  he  generated  that  are  orthogonal  over 
a  finite  set  of  discrete  values  of  the  variable  x.  Note  that  these  values 
need  not  be  equally  spaced^  a  condition  that  is  obviously  necessary  for  stress- 
rupture  data. 


Scaling  of  Polynomial  Argment 

From  the  recurrence  relation  (A3)  vith  Qj  =  1^  it  follows  that  the  lead¬ 
ing  term  of  is  xK  Therefore^  depending  on  the  values  of  Xj_^  the 

values  of  become  very  large  or  veiy  small.  This  procedure  can 

lead  to  a  loss  of  significant  figures  in  performing  the  calculations.  It  is 
shown  in  reference  7^  by  comparison  with  the  Chebyshov  polynomials^  that  if 
X  is  scaled  so  that  all  the  values  of  lie  between  2  and  -2^  the  polynomial 
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values  Q.(Xj_)  vill  all  "be  of  approximately  uniform  size.  To  perform  this 
scaling^  iet  x^ax  maximum  value  of  log  a  and  minimum 

value  of  log  a;  then  let 


and  solving 

It  has 
does  indeed 


X  =  A  log  0  +  B 

(A8) 

2  =  ® 

(A9a) 

-2  =  Axjrin  +  B 

(ASh) 

for 

A  and  B  results  in  equations  (91). 

been  found  in  practice  that  scaling  the  values  of  x 
preserve  the  significance  of  the  calculations. 

as  indicated 

Least-Squares  Procedure 

In  terms  of  the  orthogonal  polynomials^  equation  (s)  can  he  written 


m+1 

y  =  ^ 

0=1 

To  find  the  hest  values  of  y^^  and  Uj  that  fit  the  data^  the  sum  of  the 
squares  of  the  residuals  is  minimized.  Thus  let 


n 


j=i 


^0*^0 


(Xi) 


(All) 


Then  in  order  to  find  the  values  of  yg^  and  uj  that  will  make  S  a  minimum^ 

S  is  differentiated  in  turn  with  respect  to  and  each  u^^  and  the  resulting 

equations  are  set  equal  to  zero.  When  this  is  done^  the  following  set  of  equa¬ 
tions  is  obtained: 


^^a  ^1^1  ^2^2  •  • 

a^Ya  +  +  0  +  . 

EgYa  +  0  +  b2U2  +  . 


+  ^+1%+1  =  cq 

•  °  =  *^1 
.  +  0  =  C2 


(A12) 


m 


+iya 


+  0  +  0  + 


\+l%+l  “  °m+l 


where 
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=0  =  t 

i=l 

n 

^ ^  J  -  1^2  .  .  .  m  +  1 

i=l 
n 

bj  =  ^  T?Qj(Xi)  j  =  1,2  .  .  .  m  +  1 

i=l 

°0  =  S 

1=1 

n 

c.  =  >  T.y.Q.(X.  )  j  s=  1^2  .  .  .  m  +  1 

J  /  j  1  1  J  1 

1=1 


r  (^3) 


It  Is  to  be  noted  that  the  only  nonzero  elements  In  the  coefficient  matrix  of 
equations  (A12)  are  the  diagonal  elements  and  the  elements  of  the  first  row  and 
first  column.  All  the  other  elements  are  zero  because  of  the  orthogonality 
properties  of  the  polynomials  used.  This  Is  one  of  the  major  advantages  In 
using  orthogonal  polynomials.  In  the  usual  case  of  data  fitting^  all  the  ele¬ 
ments  of  the  first  row  and  first  column^  except  for  the  first  element^  would 
also  be  zero;  and  the  equations  would  be  completely  uncoupled^  each  uj 
being  computed  completely  Independent  of  the  others^  without  the  necessity  of 
solving  any  sets  of  equations  with  the  resultant  loss  of  significant  figures. 

In  this  particular  case  because  of  the  added  constant  y^^  the  equations  are 
not  completely  uncoupled^  but  they  are  very  nearly  uncoupled  and  can  readily 
he  solved.  Thus  for  any  equation  after  the  first 


Substituting  into  the  first  equation  and  solving  for  ya  give  immediately 

m+1 

^0 


(A14) 


ya  = 


'  -H 

J=1  ^ 


m+1  2 

j=i  J 


(A15) 
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1  68 

PRMT 

169 

PRMT 

17C 

PRMT 

171 

PRMT 

1  72 

PRMT 

173 

PRMT 

174 

PRMT 

175 

PRMT 

176 

PRMT 

17/ 

PRMT 

17B 

PRMT 

I  79 

PRMT 

IbO 

PRMT 

lol 

PRMT 

IBz 

PKMT 

183 

PRMT 

lb4 

PKMT 

185 

PR^:T 

1  86 

PRMT 

16  7 

PRMT 

iba 

PRMT 

189 

pkM  r 

190 
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DO  300  K3=l,KrA 

c 

c  CALCULATE  TAU  ANu  rAU»»2 

C 

DO  120  K=1,N 

rjIFF=ABS( TT(K)-TTAIK3) ) 

If  (TUIFF)  116,115,118 
115  rAU{K)=0, 

CO  rO  119 

118  TAU(K) =SIGQ(K)*TaiFF»*R(K4) 

119  rAUSQR(K)  =  rAU{K)»»2 

120  lQNTIISiuE 

c 

C  EVALUATE  PULY.mOMIALS 

c 

CALL  PaLY(MAX) 

C 

UQ  200  K2=1,KM 

C 

C  UETERMINL  Y(A) 

C 

CALL  YSUBA  (l^(K2)  ) 

C 

C  CALCULATE  THEORETICAL  LOG  TIMES  ANO  TIMES 

C 

CALL  YTH(M(K2 ) ) 

C 

C  COMPUTE  RESIDUAL 

C 

CALL  R£SID(M(K2) ) 

C 

C  MAKE  OiNiE  ENTRY  IN  SUMMARY  TABLE 

C 

NTRY=NTRY+1 
TABLE! i,NTRY)=0(K5) 

TABLE(2,NTRY)=R(K4) 

ITBLE! 3,NTKY)=M(K2) 

TABLE! 4,NTRY)=TTA(K3) 

TABLE!  5,'NTRY)=^YA 
TABLE! 6,NTRY)=00 
rRGGR2=NTRY.E0.2*L INES 

IF  ! TR0GR2)  GO  TO  170 
CO  TO  190 
C 

C  OUTPUTS  ONE  PAGE  OF  SUMMARY  TABLE 

C 

c 

C  OUTPUT  TITLE  AND  HEADINGS  FOR  SUMMARY  TABLE 

C 

170  WRITE  16,9005)  ( T I TLE ( K ) , K= 1 , I  2 ) , DEGREE 

IF  (LINES. EQ,51)  WRITE  (6,9006)  KTA , KR , KO , KM , N 
WRITE  (6,9007) 

TRCCR1=NTRY.LE.L INES 
LIMIT=LINES 

IF  (TRCGRi)  LIMIT^NTRY 
uQ  180  K=l, LIMIT 

WRITE  (6,9008)  (TABLE! I, K), 1=1, 2), I TBLE( 3,K) , ( TABLE! I, K) , 1=4,6) 
IF  ( TRCGRI )  GO  TO  180 
K0L2=K+L INES 

IF  ( TRGGR2)  GO  TO  175 

IF  ( K0L2.GT.i\lTRY)  GO  TO  180 

175  WRITE  (6,9009)  ( TABLE !  I , KOL 2 ) ,  I  =  1 , 2 ) , I TBL E ( 3 , KUL 2 ) , 

I  ( TABLE!  I,K0L2) , 1=4,6) 

180  CONTINUE 

NTRY=0 
L  INES  =  55 


PRMT  191 
PRMT  192 
PRMT  193 
PRMT  194 
PRMT  195 
PRMT  196 
PRMT  197 
PRMT  198 
PRMT  199 
PRMT  200 
PRMT  201 
PRMT  202 
PRMT  203 
PRMT  204 
PRMT  205 
PRMT  2u6 
PRMT  207 
PRMT  208 
PRMT  209 
PRMT  210 
PRMT  211 
PRMT  212 
PRMT  2i.3 
PRMT  214 
PRMT  215 
PRMT  216 
PRMT  217 
PRMT  218 
PRMT  219 
PRMT  220 
PRMT  221 
PRMT  222 
PRMT  223 
PRMT  224 
PRMT  225 
PRMT  226 
PRMT  227 
PRMT  228 
PRMT  229 
PRMT  230 
PRMT  231 
PRMT  232 
PRMT  233 
PRMT  234 
PRMT  235 


PRMT  236 
PRMT  237 
PRMT  238 
PRMT  239 
PRMT  240 
PRMT  241 
PRMT  242 
PRMT  243 
PRMT  244 
PRMT  245 
PRMT  246 
PRMT  24/ 
PRMT  248 
PRMT  249 
PRMT  250 
PRMT  251 
PRMT  252 
PRMT  253 
PRMT  254 
PRMT  255 
PRMT  256 
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IF  (TftGGR3)  GO  TO  1000 
C 

C  SAVE  VALUES  PRODUCING  SMALLEST  RESIUUAL 

C 

190  IF  (DDl.LE.OO)  GO  TO  200 

Ml  =  M(K2) 

TTAl  =  TTAU3) 

R1  =  R(K4) 

Ql  =  QIK5) 

TAl  =  YA 
DD1=D0 

200  CONTINUE 

300  CONTINUE 

400  CONTINUE 

500  CONTINUE 

TRGGR3=.TRUE. 

IF  (NTRY.NE.O)  GO  TO  170 
C 

c  END  MAJOR  LOOP 

C 

C  OUTPUT  OPTIMUM  VALUES  AND  HEADING  FOR  FULL  TABLE 

C 

iOOO  CONTINUE 

1010  vJRlTL  (6,90CS)  (  T  I  TL  E  I  K  )  ,  K=  1  ♦  12  )  i  DEGREE 

LINFS=3 

1020  WRITE  (6,9010)  Q 1 , R  1 , M 1 , T T A  I , YA I , DD I 
LINES=LINES+S 

1030  WRITE  (6,9011) 

LlNES=LlNES+3 

C 

c  CALCULATE  THEORETICAL  TIMES,  RATIOS  UF  DIFFERENCE^ 

C  TO  RESIDUAL,  AND  VALUES  OF  THE  PARAMETER,  FOR  THE 

C  PARAMETER  PRODUCING  THE  MINIMUM  RESIDUAL 

C 

DO  1035  K=l,N 
T0IFF  =  ABS(TT(K)-'TTAI) 

SIGQ(K)=SIGMA(K)**Q1 

IF  (TDIFF)  1032,1031,1032 

1031  TAU(K)=0. 

UO  TO  1034 

1032  TAU{K)=SIGG(K)»TDIFF*»R1 

1034  TAUS(dR(K)  =  TAU(K)*»2 

1035  CONTINUE 
DD=DDl 

CALL  PULY(Ml) 

CALL  YSUBA(Ml) 

CALL  YTh  (Ml) 

CALL  RATIOl 
CALL  PARAM 
C 

C  OUTPUT  FULL  TABLE 

C 

K  =  0 

1040  K  =  K+1 

WRITE  (6,9012)  TT ( K ) , S I GMA ( K ) , X ( K ) , T ( K ) , T I  ME ( K ) , Y ( K ) , YY ( K ) , 
1  RATIOl K) ,PP(K) 

L1NES=L INES+1 

IF  (K.EQ*N)  GO  TO  1 
IF  (LInES.lt. 60)  GO  TO  1040 
WRITE  (6,9005)  ( T  ITLE ( KKK ) , KKK= 1 , 12 ) , DEGREE 
WRITE  (6,9011) 

LINES=6 

GO  TO  1040 
C 

C  END  OF  PROGRAM 

C 


PRMT  2l>7 
PRMT  258 
PRMT  259 
PRMT  260 
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PRMl  262 
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PRMT  264 
PRMT  265 
PRMT  266 
PRMT  267 
PRMT  268 
PRMT  269 
PRMT  270 
PRMT  2/1 
PRMT  2/2 
PRMT  2/3 
PRMT  2/4 
PRMT  275 
PRMT  276 
PRMT  277 
PRMT  278 
PRMT  279 
PRMT  280 
PRMT  281 
PRMT  282 
PRMT  283 
PRMT  284 
PRMT  285 
PRMT  286 
PRMT  28/ 
PRMT  288 
PRMT  289 
PRMT  290 
PRMT  291 
PRMT  292 
PRMT  293 
PRMT  294 
PRMT  295 
PRMT  296 
PRMT  297 
PRMT  298 
PRMT  299 
PRMT  300 
PRMT  301 
PRMT  3C2 
PRMT  T03 
PRMT  304 
PRMT  305 
PRMT  306 
PRMT  307 
PRMT  3u8 
PRMT  309 
PRMT  310 
PRMT  311 
PRMT  3i2 
PRMT  313 
PRMT  314 
PRMT  315 
PRMT  316 
PRMT  317 
PRMT  318 
PRMT  319 
PRMT  320 
PRMT  321 
PRMT  322 
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C  FORMAT  STATFMENTS  FOR  PROGRAM 

L 

C  FURMATi>  FOR  IMPUT 

C 

9001  FORMAT  (12A6) 

9002  FORMAT  (A3,F10.0) 

9003  FORMAT  (A3, 12) 

9004  FORMAT  ( A3 , OPF 1 0 . 0 , 3P F 1 0 . 0 , CPF  1 0 . 0 ) 

:  FORMAT:>  FOR  OUTPUT 

:  TITLE  (SKIPS  TO  NEW  PAGE) 

9003  FORMAT! 1H1,2UX, 12A6/1H  ,30X,A6,10H  PARAMETER/lH  ) 

:  SUMMARY  OF  INPUT 

9006  FORMAT  ( IH  , lOX , 4 5HCREEP/RUPTURE  PARAMETERS  ARF  INVESTIGATED  FOR/ 
IIH  ,I2,18H  VALUE(S)  OF  T(A),,I3,25H  TEMPERATURE  E XPUNEN T ( S ) , , 1 3 , 
224H  STRESS  EXPONENTIS),  AND,I3,14H  POL YNOM I AL ( S ) / 1  f  1  , lOX , 5HUS I NG, 
314, 12H  DATA  POINTS/IH  ) 

:  HEADINGS  FOR  SUMMARY  TABLE,  ONE  LINE  OF  SUMMARY  TABLE 

9007  FORMAT  ( IH  , 2 ( 2 X , 1 Hg , 7 X , 1 HR , 6X , IHM , 5X , 4HT ( A ) , 5 X , 4 H Y { A ) , 4X , 

1  8HSTU.DEV. , i0X)/lH  ) 

9008  FORMAT  ( IH  , OPF 5 . 2 , F8 . 2 , 1 5 , F9 . 0 , F 1 0 • 2 , IPE 11 . 2 ) 

9009  FORMAT  (lH+,S8X,0PF5.2,Fe.2,I5,F9.0,F10,2,lPE11.2) 

OPTIMUM  VALUE:> 

9010  FURMAKIH  10X44HVALUES  PRODUCING  SMALLEST  STANDARD  DE V  I A T 1  ON/ 3H0g^ 
1F3*2,4H,  R  =  F!>,2,4H,  M=I2,7H,  T  {  A  )  =  F6 . 0 , 7H  ,  Y  {  A  )  =F  9  •  3 , 1 1 H  ,  STD. DEV 
2=1PE9.2/1H0) 

HEADINGS  FOR  FULL  TABLE,  ONE  LINt  OF  FULL  TABLE 

9011  FORMAT  ( 5H  TEMP , 4X , 6M S TRE S S , 3 X , 3HL0G, 6X , 4HT I ME , 3 X , 6HC AL C 1 D , 3 X , 
13HL0G, 3X,8HCALC  LOG , 2 X6HUE V/ Su , 3X , 9HPARAMETER/ IH  , 8 X , 6 H ( * E -3 ) , 2 X , 
26H$TRESS, i4X,4HTIMe, 5X,4hTIME,4X,4HT1ME/IH  ) 

9012  FORMAT  ( IH  , OPF 5 . 0 , F8 . 1 , F8 . 3 , 2 FI 0 .  1 , 3F9 . 3 ,  1  P F  I  2 . 3 ) 

9999  FORMAT  ( IHl ) 
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^OOO  o  u  u  u  ^oo 


IBFIC  PRMtiLK  LIST, REF, DECK 

SETS  FIRST  PULYNOMIAL  TO  UNITY  AT  ALL  STATIONS  AND  STORES  PRMB  1 

ALPHAMERIC  CODE  WORDS  PRMB  2 

PRMB  3 

BLOCK  DATA  PRMB  4 

COMMON  /PLYNML/QQ( 21, 200 ) , OTHERS! 85 ) /END/ END/ NAME S/NAMES ( 2 )  PRMB  5 

DATA  (QQ{1,K),K=1,200)/200»1./,END/3HEND/,  PRMB  6 

1  (NAMES{K),K=1,2)/12HLAR$0NLIN£AR/  PRMB  7 

EMO  PRMB  8 


IBFTC  PARAM  L I  ST, REF , DECK 

SUBROUTINE  FOR  EVALUATING  THE  PARAMETER  AT  EACH  POINT  PARM  I 

PARM  2 

SUBROUTINE  PARAM  PARM  3 

PARM  4 

COMMON  /FDATA/SIGa(200),TAU(200),OTHERS{600),Y(200)  PARM  5 

1  /CALC/PP(200),OTHERl(600)/N/N  PARM  6 

2  /PLYNML/0THER2(4221) ,YA,0THER3(63)  PARM  7 

PARM  8 

DO  10  K=1,N  PARM  9 

PP(K)  =  (Y(KJ-SIGQ(K)*YA)/TAU{K)  PARM  10 

10  CONTINUE  PARM  11 

return  PARM  12 

end  PARM  13 


IBFTC  YTH  LIST, REF, DECK 

SUBROUTINE  FOR  CALCULATING  TIMES  AND  LOG  TIMES  FROM  THE  PARAMETER  YTH 


YTH 

SUBROUTINE  YTH(M)  YTH 

C  YTH 

COMMON  /CALC/OTHERS(400) ,TIME(200) , YY{200)  YTH 

1  /FDATA/SIGQ! 200) , TAU ( 200 ), OTHER  1 ( 800 )  YTH 

2  /PLYNML/QQ(21,200) ,U(21) ,YA,QTHER2{ 63)  YTH 

3  /N/N  YTH 

C  YTH 

DO  10  K=1,N  YTH 

YY{K)  =  0,  YTH 

10  CONTINUE  YTH 

Ml  =  M+1  YTH 

DO  30  K=1,N  YTH 

DO  20  J=1,M1  YTH 

YY(K)  =  YY(K)+QQ( J,K)*U( J)  YTH 

20  CONTINUE  YTH 

YY(K)  =  TAU(K)*YY(K)+SIGQ(K)*YA  YTH 

TIME(K)  =  10.*»YY(K)  YTH 

30  CONTINUE  YTH 

RETURN  YTH 

END  YTH 


10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 


23 


oonoorooooono^y*  n  o  oon<^ 


IBFTC  RATIOl  L I  ST , REF, DECK 

SUBROUTINE  FOR  CALCULATING  RATIOS 

OF  INDIVIDUAL  RESIDUALS  TO  ROOT-MEAN-SQUARE  RESIDUAL 
SUBROUTINE  RATIOl 

COMMON  /FDATA/OTHERSI 1000) ,Y( 200) 

1  /CALC/OTHERK 200) , RAT  10(200 ) ,OTHER2( 200) ,YY( 200) 

2  /N/N/DD/OD 

DO  10  K=1,N 

RATIO(K)  =  ABS(Y(K)-YY(K))/DD 
10  CONTINUE 

RETURN 

END 


RATO  1 
RATO  2 
RATO  3 
RATO  4 
RATO  5 
RATO  6 
RATO  7 
RATO  8 
RATO  9 
RATO  10 
RATO  11 
RATO  12 
RATO  13 
RATO  14 


IBFTC  RESIO  LIST, REF, DECK 

SUBROUTINE  FOR  CALCULATING  RESIDUAL 

RESD 

THE  RESIDUAL  IS  BASED  ON  THE  LOG  OF  THE  TIME.  RESD 

IT  IS  DEFINED  AS  THE  SQUARE  ROOT  OF  THE  SUM  OF  THE  SQUARES  OF  RESD 

THE  INDIVIDUAL  RESIDUALS  DIVIDED  BY  THE  DIFFERENCE  BETWEEN  THE  NUMRESD 
BER  OF  DATA  POINTS  AND  THE  DEGREES  OF  FREEDOM.  THE  DEGREES  OF  RESD 
FREEDOM,  KK,  DEPENDS  ON  THE  PARAMETER  (SEE  MAIN  BODY  OF  REPORT).  RESD 

KK=2  FOR  LARSGN-MILLER  PARAMETER  RESD 

KK=3  FUR  LINEAR  PARAMETER  RESD 

KK=5  FOR  GENERAL  PARAMETER  RESD 

RESD 

DD  =  SURTI ( Y-YY)**2/( N-M-KK) )  RESD 

RESD 

SUBROUTINE  R£SID(M)  RESD 

Q  RESD 

COMMON  /FDATA/OTHERSI 1000) , Y{200)  RESD 

1  /CALC/OTHER1(600) ,YY( 200)  RESD 

2  /DD/DD/N/N/DEGREE/DEGREE/NAMES/FAMES(2)  RESD 

C  RESD 

IF  (DEGREE. EQ.FAMES(2) )  GO  TO  20  RESD 

IF  (DEGREE. EQ. FAMES( 1) )  GO  TO  10  RESD 

KK  =  5  RESD 

GO  TO  30  f^ESD 

10  KK  =  2 

GO  TO  30  RESD 

20  KK  =  3  <^ESD 

30  D  =  N-M-KK  RtSO 

DO  =  0.  RESD 

DO  40  K=1,N  RESD 

OD  =  D0+( Y(K)-YY(K) )**2  RESD 

40  CONTINUE  J^ESO 

DO  =  SQRT(DO/D)  ^^ESD 

RETURN  RESD 

end  RESD 
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$IBFTC  YSUBA  L I  ST , RFF t DECK 


c 

SUBROUIINE  FUR  EVALUATING  Y(A) 

YSUB 

1 

c 

YSUB 

2 

c 

THIS  Subroutine  also  evaluates  the  QUAiMTiTiESt  Ut  necessaky 

YSUB 

3 

c 

FOR  DETERMINING  THE  THEORETICAL  LOG  TIMES. 

YSUB 

4 

c 

YSUB 

5 

SUBROUTINE  YSUBA(M) 

YSUB 

6 

c 

YSUB 

/ 

COMMON  /PLYNML/UQ{ 2If 200) tU( 21 ) ,YA,A(21),B(21)»C(2l) 

YSUB 

8 

1  /FOATA/SIGQ( 200) ,TAU( 200) , TAUS0R( 200) , OT HER S ( 400 ) , 

YSUB 

9 

2  Y(200)/N/N 

YSUB 

10 

c 

YSUB 

11 

AO  =  0. 

YSUB 

12 

CO  =  0. 

YSUB 

13 

DO  10  K=L,N 

YSUB 

14 

AO  =  AO+SIGg(K) **2 

YSUB 

15 

CO  =  CO+SIGQ(K) *Y(K) 

YSUB 

16 

10 

CONTINUE 

YSUB 

17 

Ml  =  M+1 

YSUB 

18 

DO  20  J=l,Ml 

YSUB 

19 

A(J)  =  0. 

YSUB 

20 

B{ J)  =  0. 

YSUB 

21 

C{J)  =  0. 

YSUB 

22 

20 

CONTINUE 

YSUB 

23 

DO  40  J=l,Ml 

YSUB 

24 

DO  30  K=1,N 

YSUB 

25 

AIJ)  =  A(  J)^':>IGg(K)*TAU(K)*gOUtK) 

YSUB 

^6 

6(J)  =  8(  J)^-TAUSQR(K)*Qg(  J,K)*»2 

YSUB 

27 

CIJ)  =  C(J)  +  TAU{K)»Y{K)*gg( J,K) 

YSUB 

28 

30 

CONTINUE 

YSUB 

2  9 

40 

CONT INUE 

YSUB 

30 

SUMl  =  0. 

YSUB 

31 

SUM2  =  0. 

YSUB 

32 

uO  50  J=ItMl 

YSUB 

33 

AOB  =  A(J)/B(J) 

YSUB 

34 

SUMl  =  SUMH-AOB*C{  J) 

YSUB 

35 

SUM2  =  SUM2  +  A0B»A( J  ) 

YSUB 

36 

50 

CONTINUE 

YSUB 

37 

YA  =(CO“SUMl)/( A0'$UM2) 

YSUB 

38 

DO  60  J=1,M1 

YSUB 

39 

U{J)  =  (C{  J)*-AI  J)*YA)/i3(  J) 

YSUB 

40 

60 

CONTINUE 

YSUB 

41 

RETURN 

YSUB 

42 

LND 

YSUB 

43 

25 


$lbFTC  PULY  LISTtftEFtOhCK 


c 

SUBROUTINE  FOR  EVALUATING  ORTHOGONAL  POLYNOMIALS 

POLY 

1 

c 

POLY 

2 

c 

ALL  POLYNOMIALS  UP  TO  MAXIMUM  DESIRED  DEGREE  ARE  EVALUATED 

POLY 

3 

c 

AT  EACH  DATA  POINT 

POLY 

4 

c 

POLY 

5 

c 

THE  FIRST  POLYNOMIAL  IS  IDENTICALLY  EQUAL  TO  UNITY 

POLY 

6 

c 

THESE  VALUES  ARE  STORED  BY  A  BLOCK  DATA  SUBROUTINE 

POLY 

7 

c 

POLY 

8 

SUBROUTINE  POLYIM) 

POLY 

c 

POLY 

10 

COMMON  /FDATa/OTHERU 400) ,TAUSQR(200) ,0THER2(200) ,XX{200) , 

POLY 

il 

1  OTHER3(200) 

POLY 

12 

2  /PLYNML/QQ( 21,200) , OTHERS (45) , ALPHA ( 20 ), BE T A ( 20 ) 

POLY 

13 

3  /N/N 

POLY 

14 

c 

POLY 

15 

Si  =  0. 

POLY 

16 

S2  =  0- 

POLY 

17 

DO  10  K=1,N 

POLY 

18 

Si  =  Si+XX(K)frTAUSQR(K) 

POLY 

19 

S2  =  S2+TAUSuR(K) 

POLY 

^0 

10 

CONTINUE 

POLY 

21 

ALPHAd)  =  SI/S2 

POLY 

22 

DO  20  K=1,N 

POLY 

23 

QQ(2,K)  =  XX(K)“ALPHA( 1) 

POLY 

24 

20 

CONTINUE 

POLY 

25 

IF  (M.LE.l)  RETURN 

POLY 

26 

DO  50  K=2,M 

PULY 

27 

SI  =  0, 

POLY 

28 

S2  =  0. 

POLY 

29 

S3  =  0* 

POLY 

30 

S4  =  0. 

POLY 

31 

DO  30  J=1,N 

POLY 

32 

D1  =  TaUSQRI J)*QW{K, J) 

POLY 

33 

02  =  Dl*gO(K,J) 

POLY 

34 

SI  =  Si+XX(J)*D2 

POLY 

35 

S2  =  S2+D2 

POLY 

36 

S3  =  S3  +  XX(  J)*Dl*Qg(K‘-l,  J) 

POLY 

37 

S4  =  S4+TAUSUR( J)*Q0(K“1,J)*»2 

POLY 

38 

30 

CONTINUE 

PULY 

39 

ALPHA(K)  =  SI/S2 

POLY 

40 

BETA(K)  =  S3/S4 

POLY 

41 

DO  40  J=1,N 

POLY 

42 

QQ(K+1,J)  =  (XX(J)-ALPHA(K))*QQ(K,J)-BETA(K)*QQ{K-1,J) 

POLY 

43 

40 

CONTINUE 

POLY 

44 

50 

CONTINUE 

POLY 

45 

RETURN 

POLY 

46 

END 

POLY 

47 

n  o  o  o  ^ 


IBFTC 

C 

C 

10 

20 


SCALE  L IST,REF,0ECK 
SUBROUTINE  FUR  SCALING  LOGS  OF  STRESS 

THt  SCALED  VALUES  LIE  IN  THE  REGION  -2  TO  2 

SUBROUTINE  SlALE 

COMMON  /FDATA/QTHERU 600) ,X(200) » XX (200) , OTHER? ( 200 ) /N/N 

BIG  =  0. 

SMALL  =  U£5 

UO  10  K=ltN 
BIG  =  AMAXl(bIG,X(K)  ) 

SMALL  =  AMIN1(SMALL,X(K) ) 

CONTINUE 

A  =  4W(BIG-SMALL) 

B=2.*(blG+SMALL)/(BIG-SMALL) 

UO  20  K=l,N 
XXIK)  =  A»X(K)-6 
CONTINUE 
RETURN 
END 


SCAL  i 
SCAL  2 
SCAL  3 
SCAL  4 
SCAL  S 
SCAL  6 
SCAL  7 
SCAL  a 
SCAL  9 
SCAL  10 
SCAL  11 
SCAL  12 
SCAL  13 
SCAL  14 
SCAL  15 
SCAL  16 
SCAL  17 
SCAL  18 
SCAL  19 
SCAL  20 
SCAL  21 
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TABLE  I.  -  CALCULATION  OF  POLYNOMIALS  FOR  THEORETICAL  DATA  FOR  THIRD  DEGREE  POLYNOMIAL 


[Temperature  intercept^  T^,  600°  F.  ] 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

Index, 

i 

Tempera¬ 

ture, 

Oji 

Time, 

hr 

Stress, 

psi 

log  t 

log  a 

a^(T-Taf 

Scaled 
log  a, 

X 

Polynomial 

% 

Q4 

1 

1100 

4954. 68 

56  300 

3.  69501 

4.  75051 

500 

2.0 

1 

1. 3619 

-0. 19594 

-0. 50845 

2 

1100 

11365. 9 

19  800 

4. 05560 

4.  29666 

500 

1.  3806 

1 

.  30341 

-1. 6875 

-.57205 

3 

1200 

625. 342 

30  300 

2. 79612 

4.  48144 

600 

1. 6328 

1 

. 85576 

-1. 4583 

1.  7195 

4 

1200 

2908. 

5  080 

3. 46359 

3.  70586 

600 

.57433 

1 

-.80869 

1. 5741 

-1.  2980 

5 

1300 

117. 371 

12  900 

2.06956 

4. 11059 

700 

1.1267 

1 

. 18925 

2.5067 

-1.0745 

6 

1300 

1340. 

778 

3. 12710 

2.  89098 

700 

-.53777 

1 

-2.2709 

-.  90880 

.  92564 

7 

1400 

34.  4856 

4  190 

1.53764 

3.  62221 

800 

. 46017 

1 

2. 8030 

.  015445 

-.77354 

8 

1400 

995.25 

66 

2. 99793 

1.  81954 

800 

-2.0 

1 

.51879 

78251 

98133 

TABLE  II.  -  INTERMEDIATE  CALCULATEONS  FOR  THEORETICAL 
DATA  FOR  THIRD  DEGREE  POLYNOMIAL 
[Temperature  intercept,  Tg^,  600°  F.  ] 


1 

2 

3 

4 

5 

6 

7 

Index, 

j 

a 

f3 

a 

b 

c 

u 

0 

8.0 

,23.  743 

1 

0.27092 

0. 

5200. 

3. 48  XlO® 

14897. 

-9. 9146X10“^ 

2 

-.57813 

1. 6548 

786. 18 

5.  7589 

2616. 

-8. 4266X10“^ 

3 

.  28432 

1.  3315 

465. 68 

7.6678 

3796.  9 

-8. 1771X10“^ 

4 

. 41260 

. 80618 

286.  01 

6. 1816 

2694. 6 

-3. 6513X10-6 
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TABLE  III.  -  PIT  FOR  SEVERAL  VALUES  OF  LINEAR 


PARAMETER  FOR  THEORETICAL  DATA 


Degree  of 
polynomial 

Temperature^ 

Ta 

Variable, 

ya 

Deviation 

2 

500 

10.54 

0. 008049 

3 

500 

10.  54 

. 009786 

4 

500 

10.55 

. 010660 

2 

600 

9.49 

.  004859 

3 

600 

9.50 

.  000002 

4 

600 

9.50 

.  000003 

2 

700 

8.  44 

.  015412 

3 

700 

8.  46 

. 013937 

4 

700 

8.45 

.  014469 

TABLE  IV.  -  COMPOSITION  OF  STEELS  RECEIVED 
FROM  GERMAN  COOPERATIVE  LONG¬ 
TIME  CREEP  PROGRAM 

[As -received^  SO-mra-diam.  bar  stock.  ] 


Element 

Composition^  percent 

Steel 

c 

(23b  CK) 

p 

(I4a  PA) 

K 

(27b  KK) 

Carbon 

0.065 

0.270 

0.068 

Silicon 

.47 

.26 

.45 

Manganese 

.60 

.60 

.73 

Chromium 

17.24  ^ 

2.62 

16.14 

Molybdenum 

2.08 

.27 

2.10 

Columbium  and 

.02 

Trace 

.44 

tantalum 

Nickel 

11.  90 

.14 

13. 12 

Titanium 

.39 

Trace 

Trace 

Vanadium 

.10 

.26 

.05 

Tungsten 

Less  than 

Trace 

Trace 

0.005 

TABLE  V.  -  NASA  RUPTUEE  DATA 


(a)  Steel  K  (27b  KK) 


Temperature  ^ 
T. 

Op 

Stress^ 

o, 

psi 

Tirae^ 

hr 

Temperature^ 

T. 

Op 

Stress^ 

psi 

Time^ 

hr 

Temperature^ 

T. 

Op 

Stress^ 

0, 

psi 

Time; 

hr 

1022. 00 

77  000.000 

1.500 

1600. 00 

20  000.000 

0.400 

^1112.  00 

60  000.000 

12.  900 

^1022.00 

72  500.000 

13.  800 

1560. 00 

20  000.000 

1.900 

2-1110.  00 

60  000.000 

34. 

^1022. 00 

72  000.000 

10. 

1520.00 

20  000.000 

4.450 

2-1080.  00 

60  000.000 

52.  200 

^1022.00 

70  000.000 

36.  700 

^1480,  00 

20  000.000 

23. 700 

2-1080.  00 

60  000.000 

37.  400 

^-1022.  00 

68  000.000 

60. 400 

^1460.  00 

20  000,000 

25.500 

3-1050. 00 

60  000.000 

239. 

^■1022. 00 

66  000.000 

73.  300 

®1440.00 

20  000.000 

38. 

ai030.00 

60  000.000 

445. 

^1022. 00 

66  000.000 

107.  600 

^400. 00 

20  000.000 

136. 800 

2-1022.  00 

60  000.000 

989.  900 

^^102 2.  00 

65  000.000 

201. 300 

^1360. 00 

20  000.000 

394. 800 

^10 20. 00 

60  000.000 

817.500 

^•102 2. 00 

62  500.000 

250. 400 

^1340. 00 

20  000.000 

704.  600 

1040. 00 

75  000.000 

.330 

^1022.00 

60  000.000 

990. 

^1320.  00 

20  000.000 

1212. 

1022.00 

75  000.000 

5.850 

®1022.  00 

60  000.000 

817.500 

1320.00 

40  000.000 

2.700 

^-1000. 00 

75  000.000 

15.600 

^102 2. 00 

55  000.000 

3  680. 

1290.00 

40  000,000 

7.500 

2-980. 00 

75  000.000 

46.  500 

1112. 00 

68  000.000 

.  750 

®1260.00 

40  000,000 

15. 200 

2-960.  00 

75  000.000 

138. 

1112. 00 

65  000.000 

2.250 

^1230.00 

40  000.000 

44.  400 

^940. 00 

75  000.000 

542. 

1112.  00 

62  500.000 

4.300 

2-1170.  00 

40  000.000 

377. 

^920.00 

75  000.000 

579.  600 

^1112.00 

60  000.000 

13. 900 

2-1140.  00 

40  000.000 

1417. 

2-1120.00 

50  000.000 

186. 100 

®-1112.00 

57  500.000 

22.  700 

^•1125.  00 

40  000.000 

2110. 

^1200.  00 

40  000.000 

130. 200 

^1112.00 

55  000.000 

51. 500 

^1112. 00 

40  000.000 

5367. 

2-1280. 00 

30  000.000 

132. 700 

^1112.00 

52  500.000 

147.500 

1200. 00 

60  000.000 

.610 

^1340.  00 

25  000.000 

125. 800 

^1112.  00 

50  000.000 

283. 

1170. 00 

60  000.000 

1.250 

2-1500.00 

15  000.000 

51,  300 

®'1112. 00 

45  000.000 

1  020. 

1150. 00 

60  000.000 

4.  400 

2-1560.  00 

12  000.000 

41. 700 

®'1112.  00 

43  000.000 

1  579. 

1140. 00 

60  000.000 

4. 500 

21580,  00 

10  000.000 

32. 400 

1112.  00 

37  000.000 

13  140, 

2-1120. 00 

60  000.000 

10.900 

21540.  00 

10  000.000 

148. 200 

(b)  Steel  C  (23b  CK) 


Temperature, 

T, 

Op 

Stress, 

psi 

Time, 

•t, 

hr 

Temperature, 

T, 

Op 

Stress, 

psi 

tr  c+ 

^i  S 

(D 

Temperature, 

T. 

Op 

Stress, 

0, 

psi 

Time, 

hr 

21600. 00 

5  000.000 

570. 200 

21230. 00 

30  000.000 

175. 700 

^1202.  00 

36  000.000 

68.  600 

21620.  00 

5  000.000 

186.600 

21250. 00 

30  000.000 

103. 500 

^1202.00 

38  000.000 

59.  300 

21660. 00 

5  000.000 

156. 800 

21280. 00 

30  000.000 

58. 100 

21202. 00 

42  000.000 

24. 400 

^1680. 00 

5  000.000 

91. 600 

1292.00 

30  000.000 

21  300. 

21202. 00 

44  000.000 

14.500 

21700.  00 

5  000.000 

62.  700 

21310. 00 

30  000.000 

22.500 

^1202. 00 

45  000.000 

22.  900 

2i740. 00 

5  000.000 

40. 500 

2ii12.  00 

40  000.000 

667. 900 

1202.  00 

46  000.000 

7. 

21780. 00 

5  000.000 

10. 600 

21120. 00 

40  000.000 

785. 400 

1202.  00 

48  000.000 

2.850 

21425.00 

10  000.000 

1690. 

21150. 00 

40  000.000 

266. 700 

1202.  00 

49  000.000 

2.550 

21450. 00 

10  000.000 

550. 300 

21170. 00 

40  000.000 

127. 800 

1202.  00 

50  000.000 

1.470 

21480. 00 

10  000.000 

270. 

21202.00 

40  000,000 

44. 100 

^1292.  00 

18  000.000 

859.  700 

21500.00 

10  000.000 

170. 

21202.00 

40  000.000 

74. 

^1292.  00 

23  000.000 

194.  600 

21520, 00 

10  000.000 

128.  500 

21210.  00 

40  000.000 

40. 500 

^1292.  00 

25  000,000 

75, 

^1560. 00 

10  000.000 

40. 

21220. 00 

40  000.000 

37.  800 

21292.00 

28  000.000 

34. 600 

21570.00 

10  000.000 

31.500 

21240. 00 

40  000.000 

17.  200 

®-1292.  00 

29  000.000 

31. 

21600. 00 

10  000.000 

15. 800 

1270. 00 

40  000.000 

4. 500 

®-1292.  00 

32  000.000 

13. 300 

1650.  00 

10  000.000 

5.  250 

1280. 00 

40  000.000 

1.200 

21292.00 

33  000.000 

19. 800 

■  1700.00 

10  000.000 

1.  750 

1292.00 

40  000.000 

1.300 

®-1292.  00 

34  000.000 

10. 400 

21202. 00 

20  000.000 

3307. 

1300. 00 

40  000.000 

.800 

1292.00 

36  000.000 

2.750 

21260. 00 

20  000.000 

667.400 

21112. 00 

34  000.000 

2  274. 

1292.00 

37  000.000 

7.600 

21290.00 

20  000.000 

255. 

21112.00 

43  000.000 

363. 100 

1292.00 

38  000.000 

1.650 

21292.00 

20  000.000 

347. 100 

21112.00 

46  000.000 

233.  900 

^1060. 00 

60  000.000 

42.  500 

21292. 00 

20  000.000 

363. 

21112. 00 

46  000.000 

261. 400 

^1300.00 

25  000.000 

89. 600 

21320. 00 

20  000.000 

180. 400 

21112. 00 

48  000.000 

183. 100 

^1360. 00 

19  000.000 

95. 

®'1360.00 

20  000.000 

82. 

21112. 00 

50  000.000 

84.  500 

^1430. 00 

15  000.000 

71. 400 

21400. 00 

20  000.000 

28.  900 

21112.  00 

52  000.000 

65. 600 

®-1480. 00 

12  000.000 

147. 900 

1440. 00 

20  000.000 

9. 

21112. 00 

54  000.000 

39.  300 

^1570.  00 

8  000.000 

104. 

1480. 00 

20  000.000 

2.500 

21112. 00 

57  000.000 

23.  300 

^'1630.  00 

6  000.000 

140. 900 

21112. 00 

30  000.000 

4258. 

21202.00 

25  000.000 

1  074. 

^1140. 00 

34  000.000 

1077. 

®'1160,  00 

30  000.000 

1110. 

21202.00 

34  000.000 

199.  400 

^1320. 00 

15  000.000 

1505. 

21180. 00 
21202.00 

30  000.000 

30  000.000 

696. 300 
350. 

21202.00 

35  000.000 

124.  300 

^1480. 00 
^^1540.  00 

8  000.000 

6  000.000 

2237. 

1258. 

^■Data  point  used  in  parametric  analysis. 
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TABLE  V.  -  Concluded.  NASA  RUPTUBE  DATA 
(c)  Steel  P  (14a  PA) 


Temperature^ 

T, 

Opi 

Stress^ 

psi 

Time, 

hr 

Temperature, 

T, 

Op 

Stress, 

0, 

psi 

Time, 

hr 

Temperature, 

T, 

Op 

Stress, 

psi 

Time, 

hr 

932.00 

65  000.000 

3.800 

ai250. 00 

10  000.000 

19.  200 

^740. 00 

90  000.000 

57.100 

^•932. 00 

60.000.000 

14.150 

^1220. 00 

10  000.000 

42. 

a785. 00 

80  000.000 

84. 

0-932. 00 

60  000.000 

14. 400 

ailSO. 00 

10  000.000 

167. 

^•820. 00 

70  000.000 

195. 800 

O932. 00 

57  500.000 

10. 

^1170. 00 

10  000.000 

203.  400 

^■880. 00 

60  000.000 

120. 

^-932.  00 

55  000.000 

18.  900 

^"1140.  00 

10  000.000 

608. 

®'932.00 

50  000.000 

103. 500 

O932, 00 

52  500.000 

51. 

^1090. 00 

10  000.000 

2639. 

^“1022. 00 

30  000.000 

186. 700 

^-932. 00 
^*932. 00 
932.00 
1022. 00 

40  000.000 

623. 

1100. 00 

40  000. 000 

1.300 

s-1050.  00 

25  000.000 

123. 500 

30  000.000 

7  592. 

1080. 00 

40  000.000 

2.  200 

^■1090. 00 

20  000.000 

79.500 

27  000.000 

11  410. 

1060. 00 

40  000.000 

4.300 

0-1090.  00 

20  000.000 

112. 400 

58  000.000 

.580 

1050. 00 

40  000.000 

6.800 

O112O. 00 

16  000.000 

183.500 

1022. 00 

55  000.000 

.717 

1040. 00 

40  000.000 

7.400 

®'1160.  00 

13  000.000 

100. 300 

1 022. 00 

50  000.000 

1.  280 

^1020. 00 

40  000.000 

22.500 

®'1230. 00 

8  000.000 

97. 900 

1022.00 
1022.00 
1022. 00 

47  000.000 

2.450 

S-IOIO.  00 

40  000.000 

20. 100 

01290. 00 

5  000.000 

139. 700 

47  000.000 

6.  200 

aiOOO.  00 

40  000.000 

63. 300 

O74O.OO 

80  000.000 

996.600 

45  000.000 

3.500 

^990. 00 

40  000.000 

51. 200 

^*780. 00 

70  000.000 

1122. 

1022.00 

42  500.000 

6.300 

^■980.  00 

40  000.000 

80. 600 

0030. 00 

60  000.000 

948. 800 

^-1022.  00 
®'1022.00 
®'1022.00 
1415. 00 
1340. 00 
1315.00 
1290.00 

40  000.000 

22.500 

^-960. 00 

40  000.000 

192.100 

0080. 00 

50. 000. 000 

599. 

37  500.000 

12. 

^940. 00 

40  000.000 

427. 900 

0932. 00 

35  000.000 

1902. 

25  000.000 

382. 200 

^930. 00 

40  000.000 

623. 

O98O.  00 

30  000.000 

754.  800 

10  000.000 

.170 

^000. 00 

40  000.000 

2572. 

OlOOO. 00 

25  000.000 

970. 700 

10  000.000 

1. 500 

932.00 

70  000.000 

1.400 

O103O.OO 

20  000.000 

1084. 

10  000.000 

3.  700 

897.00 

70  000.000 

5.800 

O1070.00 

16  000.000 

804. 800 

10  000.000 

6.100 

^-eeo.  00 

70  000.000 

31. 200 

OllSO.  00 

01220. 00 

8  000.000 
5  000.000 

948. 500 
960. 

®Data  point  used  in  parametric  analysis. 
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TABLE  VI.  -  GERMAB  RUPTURE  DATA 


Temperature^ 

Si 

Stress j 

cr, 

psi 

Time^ 

hr 

Temperature  ^ 

Si 

°F 

Stress^ 

psi 

Time^ 

hr 

Si 

;eel  K  (27b  K 

K) 

teel  C  (23b  C 

K) 

1022.00 

76  899.999 

0.100 

1022.00 

66  899.999 

160. 

1022.00 

55  500.000 

2  000. 

1112. 00 

72  500.000 

.100 

1112.00 

64  000.000 

10. 

1112. 00 

55  500.000 

35. 

1112. 00 

44  100.000 

2  100. 

1112.00 

28  400.000 

52  000. 

steel  C  (23b  CK) 

1112.00 

14  200.000 

60  000. 

1112. 00 

17  800.000 

30  000. 

1112.  00 

28  400.000 

3  500. 

1112. 00 

28  400.000 

3  000. 

1112. 00 

28  400.000 

2  200. 

1112.00 

35  600.000 

1  200. 

1112. 00 

44  100.000 

520. 

1112.00 

51  200.000 

150. 

1112. 00 

59  800.000 

.100 

1202.00 

11  400.000 

82  790. 

1202.00 

14  200.000 

15  000. 

1202.00 

17  800.000 

6  500. 

1202.00 

22  800.000 

1  800. 

1202. 00 

28  400.000 

550. 

1202. 00 

35  600.000 

124. 

1202.00 

42  700.000 

5. 

1202.00 

52  600.000 

.100 

1292.00 

11  400.000 

30  000. 

1292.00 

11  400.000 

20  000. 

1292.00 

13  900.000 

4  500. 

1292.00 

17  800.000 

1  100. 

1292. 00 

21  400.000 

300. 

1292.00 

21  400.000 

250. 

1292.00 

27  000.000 

180. 

1292.00 

27  000.000 

140. 

1292.00 

47  000.000 

.100 

Steel  P  (I4a  PA) 


932. 00 

84  000.000 

0.100 

932.00 

75  500.000 

.100 

932. 00 

78  399.999 

2. 

932. 00 

55  500.000 

150. 

932.00 

44  100.000 

1  700. 

932. 00 

34  200.000 

2  600. 

932.00 

27  000.000 

16  000. 

932. 00 

22  800.000 

22  000. 

932. 00 

17  100.000 

100  000. 

1022.00 

72  599.999 

.100 

1022.00 

69  699.999 

.100 

1022. 00 

65  500.000 

1.  200 

1022.00 

59  800.000 

1.500 

1022.00 

35  600.000 

150. 

1022. 00 

27  000.000 

300. 

1022.00 

22  800.000 

400. 

1022.00 

22  800.000 

900. 

1022.  00 

17  100.000 

2  100. 

1022. 00 

13  900.000 

6  500. 

1022.00 

13  900.000 

8  000. 

1022.00 

11  100.000 

10  000. 

1022.00 

8  830.000 

68  000. 

TABIiE  VII.  -  CREEP  DATA  FOR  COLUMBIUM  ALLOY  FS-85 


Temperature  ^ 

'h 

°F 

Stress ; 

psi 

Time^  t^  hr 

1 -Percent  creep 

2 -Percent  creep 

5 -Percent  creep 

2005 

25  000 

0.6 

3.0 

6.1 

1900 

25  000 

26. 

33. 

45. 

1790 

25  000 

210. 

257. 

332. 

2175 

18  000 

4.  9 

7.8 

13. 

2400 

10  000 

3.4 

5.  7 

10.8 

2300 

10  000 

25.4 

41. 

68. 

2200 

10  000 

54. 

84. 

133. 

2100 

10  000 

355. 

500. 

765. 

2100 

10  000 

380 

570. 

875. 

2000 

10  000 

775. 

1325. 

2175. 

2000 

10  000 

900, 

1420. 

2000 

8  500 

2480. 

2575 

6  000 

5.6 

10. 

22.2 

2200 

6  000 

425. 

710. 

1370. 

2800 

4  000 

3,4 

6.4 

13.5 

2620 

4  000 

14.4 

26. 

56. 

2200 

4  000 

1140. 

2900 

3  000 

2.6 

5.4 

13.8 

3000 

2  000 

4.6 

9.5 

33.2 

2450 

2  000 

950. 
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”The  aeronautical  and  space  activities  of  the  United  States  shall  be 
conducted  so  as  to  contribute  ...  to  the  expansion  of  human  knowl- 
edge  of  phenomena  in  the  atmosphere  and  space.  The  Administration 
shall  provide  for  the  widest  practicable  and  appropriate  dissemination 
of  information  concerning  its  activities  and  the  results  thereof.’^ 

— ^National  Aeronautics  and  Space  Act  of  1958 
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